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ABSTRACT 

A number of modern experiments in high-energy astrophysics produce images of super- 
nova remnants (SNRs) in the TeV and GeV gamma-rays. Either relativistic electrons 
(due to the inverse-Compton scattering) or protons (due to the pion decays) may 
be responsible for this emission. In particular, the broad-band spectra of SNRs may 
be explained in both leptonic and hadronic scenarios. Another kind of observational 
data, namely, images of SNRs, is an important part of experimental information. We 
present a method to model gamma-ray images of Sedov SNRs in uniform media and 
magnetic field due to hadronic emission. These 7-rays are assumed to appear as a 
consequence of meson decays produced in inelastic collisions of accelerated protons 
with thermal protons downstream of the shock - a model would be relevant for SNRs 
without firm confirmations of the shock-cloud interaction, as e.g. SN 1006. Distribu- 
tion of surface brightness of the shell-like SNR is synthesized numerically for a number 
of configurations. An approximate analytical formula for azimuthal and radial vari- 
ation of hadronic 7-ray brightness close to the shock is derived. The properties of 
images as well as the main factors determining the surface brightness distribution are 
determined. Some conclusions which would be relevant to SN 1006 are discussed. 

Key words: ISM: supernova remnants - shock waves - ISM: pion decays - radiation 
mechanisms: non-thermal - acceleration of particles 



1 INTRODUCTION 

Cosmic rays (CRs) are an important component of the Uni- 
verse on different scales, from the Solar System to clusters 
of galaxies. Supernova remnants (SNRs), sources of galactic 
CRs, are excellent objects to study magneto- hydrodynamics 
of nonrelativistic shocks and acceleration of cosmic rays, as 
well as their mutual influence. Protons and electrons, once 
heated and accelerated on the shock, produce different types 
of emission: thermal with maximum intensity in X-rays and 
nonthermal emission in radio. X-rays and 7-rays. Experi- 
ments in high-energy astronomy observe all types of these 
emission. 

Most of galactic cosmic rays are believed to be produced 
by the forward shocks in supernova remnants (SNRs). A 
number of indirect evidences favor this expectation. In par- 
ticular, efficient proton acceleration changes the structure of 
the shock front and makes plasma more compressible that 
leads to lower adiabatic index, to increased shock compres- 
sion ratio and to some observed effects: reduced physical sep- 
aration between the forward sh ock and the "contac t discon- 
tinuity" (or reverse shock) (e.g. IWarren et aLll2005l ) or even 

J irotrusions of the ejec ta clumps beyond the forward shock 
Rakowski et al.ll201ll ): concave shape of the energy spec- 



trum (e.g. [Reynolds fc Ellisonlll992l ): growth of some turbu- 
lence modes and amplificatio n of magnetic field in the pre- 
shock region (e.g. lBeiill200i ) ; "blinking" X-ray spot s origi- 
nated from such growth of magnetic field (Uchivam a et al.l 
I2OO7I : iPatnaude fc F esen 200 71); ordered non-thermal X-ray 
strips (|Eriksen et al.i,2011i : ,Bvkov et aLlbOllI ). 

Nevertheless, there is still luck of direct observational 
confirmations that protons are accelerated in SNRs to en- 
ergies which make them responsible for 7-ray emission of 
SNRs. Broadband spectra of a number of SNRs may be 
explained either by lep tonic or hadronic scenario for 7-ray 
emission (e.g. SN 1006: lAcero et aLll201(]| ). Fermi 7-ray ob- 
servatory is expected to clarify such ambiguity, at leas t in 
bright SNRs (as in RX J1713.7-3946: lAbdo et al.ll201ll ). In 
any case, other kind of experimental information should be 
studied as well. Namely, the distribution of surface bright- 
ness in different bands caused by emission of accelerated 
particles in SNRs contains wealth of information about prop- 
erties of CRs and magnetic fields in these objects. 

A well known example is the method to estimate the 
strength of the post-shock magnetic field fro m the thickness 
of the radial profiles of X-ray brightness (e.g. lBerezhko et al.l 
I2OO3I ). Radial structure of the shock front upstream (ob- 
served in X-rays) seems to confirm back-reaction of parti- 
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cles and magnetic fi eld amplification, locally, in SN 1006 
l|Morlino et al.ll2010l ). A method to derive an aspect an- 
gle (between ambient magnetic field and the line of sight) 
in SNR from the azi muthal variat i on of t he radio bright- 
ness is developed in iPetruk et all l|2009al ) and applied to 
SN 1006 under assumption of the uniform interstellar mag- 
netic field (ISMF). More detailed consideration of the az- 
imuthal radio profile including the ISMF nonuniformity re- 
sults in constraint s on orientation of IS MF and its gradient 
around SN 1006 (|Bocchino et al.ll201ll) . Radial profiles of 
the radio brightness may const rain the time evolut ion of 
the electron injection efficiency l|Petruk et al.|[2011a^ ■ Spa- 
tially resolv ed spectral analysi s of radio and X-ray data (cis 
in SN 1006: iMicefi et al]|2009l ) determines the model, value 
and s urface variation of th e electron maximum energy in the 
SNR l|Petruk et al.ll2011al ). 

Properties of the nonthermal images of Sedov SNRs due 
to radiation of accelerated electr ons in radio. X-rays an d 7- 
rays are systematicaUy studi ed in lRevnoldd (| 19981 . [20041 ) and 
iPetruk et al] (|2009bl . Boilbl . Papers I and II respectively). 
Numerical models for synthesis of maps of adiabatic SNRs in 
uniform interstellar medium (ISM) and uniform ISMF from 
basic theoretical principles as well as approximate analyt- 
ical descriptions are developed in these papers. The main 
factors determining the azimuthal and radial variation of 
surface brightness of SNRs are determined. The role of the 
nonuniform ISM and/or nonu niform ISMF i n nontherma l 
images of SNRs are studied bv lOrlando et all (|2007l . l201lh : 
in particular, gradients of ISMF strength or ISM density 
result in several types of asymmetries in surface brightness 
distributions of SNRs in radio, hard X-ray and 7-ray bands. 
These papers are limited to the test-particle approach be- 
cause the non-linear theory of diffusive acceleration is not 
developed for shocks of different obliquity, while the obliq- 
uity dependence of various parameters is crucial in image 
modeling. 

In addition to the analysis of SNR images simulated 
from basic theoretical principles, model-independent meth- 
ods for SNR images are important. Such method for syn- 
thesis of the inverse-Compton 7-ray map of SNR from the 
radio (or hard X-ray) image and results of the spatially re- 
solved X-ra y spectral ana l ysis is developed and applied to 
SN 1006 in IPetruk erakl j2009d ). It is found that synthe- 
sized inverse-Compton 7-ray image of SN 1006 is in agree- 
ment with HESS observations. This fact favors a leptonic 
scenario for the TeV 7-ray emission of this SNR. Further 
development of this method allows us to find a new way to 
constrain the streng th of magnetic field in SN 1006 from its 
nonthermal images l|Petruk et al.ll201ld ). 

Though a leptonic scenario for 7-rays from SN 100 6 
is reasonable and promising (e.g. iPetruk et aO l2011al ). 
a hadronic origin cannot be ruled out even in view of 
the small ISM densities (upper lim it is 0.5 cm~^; e.g. 
iDubner et al.ll2002l : IXcero et al.ll200^). which are consistent 
with a hadronic scenario ( Berezhko et al.|[2009l : lAcero et al.l 
I2OIOI ). The shape of the observed spectrum in the TeV 7-ray 
range better corresp onds to the hadro nic spectrum rather 
then to leptonic one (|Acero et al.ll2010l ). In order to fit the 
observed radio. X-ray and gamma-ray emission within the 
leptonic scenario on e needs rather h i gh downstream mag - 
netic field ^ 30 /iG (|Volk et al.ll20oi : IPetruk et al.ll2011ai ): 
there is a general thought that such a field can only be pro- 



duced by efficiently accelerated cosmic ray proton compo- 
nent. 

If TeV 7-rays from SN 1006 is of hadronic origin, then 
the observed brightness map of this SNR in 7-rays with en- 
ergy > 1 TeV should reflects the distribution of protons with 
energies > 2 TeV which interact with shocked thermal pro- 
tons inside SNR. 

In the present paper, we study this possibility modeling 
7-ray images of a Sedov SNR in uniform ISM and uniform 
ISMF it would have under such scenario. Though we are 
primarily interested in comparison with SN 1006, our results 
are general and may be used for analysis of other SNRs. 



2 MODEL 

2.1 General description 

Our model closely restores model used in Papers I and II. 
Let us consider an adiabatic SNR in uniform ISM and uni- 
form ISMF . Hydrodyna mics of the remnant is given by the 
self-similar ISedovl 1 19591 ) solutions; in practice, we use their 
quite accur ate approxim ations in Lagrangian coordinates 



(Sect. 4 in Pet ruk 2000). Magnetic field is described fol- 
lowing iRevnolds ( 1998 ) ; we do not consider amplification of 
the ambient field (though its role is discussed in Sect. |4]). 
Accelerated protons are described by a number of parame- 
terizations as follows. 

At the shock, accelerated protons are distributed with 
energy as 



s+Ss{Ep) 



exp (— {Ep/Ep 

,1 



(1) 



where A'' is the differential number density of accelerated 
protons, -Bp, max is the ma^ximum energy of protons, 5s{Ep) 
reflects possible concavity of the spectrum shape as an effect 
of efficient acceleration, s and a are constants, Ks the nor- 
malization, index 's' denotes values at the shock. Different 
theoretical a pproaches results in different values of a (see 
references in iRevnolds fc Keohanelll999l : lAUen et aklboOSl : 
lOrlando et al.ll201ll : IPetruk et al.ll2011bh which are between 
1/4 and 2. The function 5s{Ep) varies a bit and very slowl y 
over decades in energy Ep (e.g. iBerezhko fc Ellisonlll999l ). 
The emission at some photon energy s is mostly determined 
by the narrow range of Ep. Therefore, 5s{Ep) produces neg- 
ligible changes to our results and we take it 5s{Ep) = 0. 
The maximum energy may vary with obliquity: 



Ep 



,it,eo) = Ep 



(t)£max(eo). 



(2) 



where -Ep,max||(i) oc V(t)', V the shock velocity, q a con- 
stant, index "||" denotes values at parallel shock, Qo is the 
obliquity angle between the ambient magnetic field and the 
shock velocity, function £max(Oo) is some function. 

There are expectations that the efficiency of injection 
(defined as the ratio of density of accelerated protons to 
density of all protons) could be most efficient at the par- 
allel shock and progressively depressed towar d regions of 
SNR surface where the shock is perpendicular ((Ellison et al.l 
1 19951 : IVolk et aLllioOsI ). In opposite, an analysis of known 
SNR images presents some hints against suc h scenario 
(|Fulbright fc Revnoldslll990l : IOrlando et al.ll200'iD : it also re- 
veals that properties of SNR images are similar in isotropic 
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(no dependence on obliquity) and quasi-perpendicular (in- 
jection prefers perpendicular shock) approaches. In our cal- 
culations, we use simple parameterizations for obliquity de- 
pendence of the injection. Namely, the obliquity dependence 
of the normalization Ks (which is proportional to the injec- 
tion efficiency) is parameterized by 



if.(t,eo) = if.|i(t)K;(eo) 



(3) 



where IC{Oo) is given by the following formulae where dif- 
ferent 'sensitivity' to the obliquity angle is given by the pa- 
rameter Qk- The quasi-parallel and quasi-perpendicular in- 
jections are represented respectively by 



ICiQo) 



exp 



exp 



(-(Go/eK) 
(-((eo-^/2)/eK) 



(4) 

(5) 



Isotropic injection IC{Oo) = const may be restored by a large 
value of 0K- Evolution of the injection efficiency may be 
accounted through K^w cc V{t)~'' where & is a parameter. 

The evolution of the proton energy spectrum in the 
SNR's interior is considered in the next subsections. 

The surface brightness is calculated integrating emis- 
sivities along the line of sight within SNR. Hadronic 7-rays 
appear as a consequence of the neutral pion and 77- meson de- 
cays produced in inelastic collisions of accelerated protons 
with thermal protons downstream of the shock; the spatial 
distribution of the target protons inside the volume of SNR 
is simply prop ortional to the local p l asma density. We use 
the approach of lAharonian fc AtovanI (|200(]| ') to calculate the 
hadronic 7-ray emissivity: 



(6) 



where and m.^- are the energy and mass of the neutral 
pion, the pion emissivity is 

CUh 



(7) 



the average energy of protons mostly responsible for creation 
of pions with energy E^r is 



Ep = mpC^ + Ett/k, 



(8) 



K — 0.17 (this value accounts also for 77-meson production), 
the cross-section of inelastic collision of a highly-energetic 
proton with a thermal proton (with negligible energy com- 
paring to the energy of the incident proton) is 



crpp(_Bp) = 28.5 -f 1.8 In {Ep/1 GeV) mb. 



(9) 



2.2 Proton energy losses due to meson production 

In order to synthesize SNR images, we need to know emis- 
sivity ((6]) in each point of the interior and, therefore, to 
describe evolution of the proton energy spectrum. Relativis- 
tic protons lose their energy due to adiabatic expansion and 
inelastic collisions. Energy losses of protons in the meson 
production is (Appendix [XJ 

{ dt 



Skcwh cTpp ( iSp ) Sp 

,kin ; 



(10) 



where c is the speed of light, nu the number density of the 
target protons, the factor 3 accounts for the production of 
tt", tt^ and 7r~ mesons. 



The losses due to proton collisions are more significant 
for larger density of the target protons and higher energy 
of incident protons. In order to feel then the proton energy 
losses are effective, let us compare them with the radiative 
energy losses of electrons 



dEc 
~dt 



rad 'J 



(11) 



where ctt is the Thomson cross-section, rrie the mass of elec- 
tron. The ratio of electron to proton losses is 



(12) 



£^p,pp "•H £'p,TcV 

where we used o-pp ~ 33 mb (|Aharonian fc AtovanI I2OO0I ) , 
_B^G magnetic field in 10~^ G, -Bo.TeV and i?p,TGV are the 
energies of electrons and protons measured in 10^^ eV. One 
can see that the losses of protons with energy 100 TeV are 
comparable to losses of electrons with energy 30 TeV in mag- 
netic field 30 ^G, if the number density of target protons is 
4 X 10* cm~'^ respectively. 



2.3 Downstream evolution of the proton energy 
distribution 

Let the energy of proton at the time ti, when it leaves the 
region of acceleration, was E^i. Then it is smaller at the 
present time t (Appendix iB]! . 



£p = £pi£ad(a)'''"'£pp(Sp,a) 



(13) 



because the terms responsible for the adiabatic fad and col- 
lisional fpp losses are equal or smaller than unity; d = a/R, 
a the Lagrangian coordinate, R the radius of SNR, 

£,i(d) = n{df/\ (14) 

where fi = n/us, index "s" denotes the value immediately 
post-shock. 



£pp{Ep,d) = (Sp/lGeV) 



l-fi(a) 



X{d), 



(15) 



^(a) and T{d) are dimensionless self-similar functions pre- 
sented in the Appendix [B] Close to the shock they behave 
as (Appendix IC)) 



£:pp(a) 



li{a) 



(16) 



where the cross section is in mb, ( = 1.21 x lO^^tsnus, 
ts = t/lOOOyrs, t the age of SNR. One can note that £pp is 
effective if the number density of the target protons or/and 
age of SNR are large. 

The energy spectrum of protons downstream of the Se- 
dov shock evolves in a self-similar way (Appendix |B} 



Np{Ep,d,t) = K{d,t)E-'' fi{d)£pp{Ep,dy 



X exp 



Ep .max £ad(a)'^'°'£pp(-E'p, a) 



(17) 



with K{d,t) = K,K{d), Kid) = a='''/2^(a)'+''(^'(^-^'''^ The 
downstream distribution of relativistic protons are modified 
mainly by the adiabatic expansion of SNR. Losses due to 
inelastic collisions affects the distribution only when C, is 
not small, that happens when tihs and/or t are large. 
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Figure 1. Hadronic 7-ray images of Sedov SNR for different models of injection: quasi-parallel (upper panel) and quasi-perpendicular 
(lower panel) for an aspect angle 0o = 90°. Parameter is 7r/2, 7r/4, 77/6, 7r/12 (from left). Hereafter, the ISMF is parallel to the 
horizontal axis and the increment in brightness is A5 = O.lSmax- 
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Figure 2. The same as on Fig. [T]for different aspect angles. Models of injection are quasi-parallel (top, Qk = ""/S), isotropic (centre) 
and quasi-perpendicular (bottom, 0/^: = 7r/6). Aspect angle <^o = 0°, 30°, 60° and 90° (from left). 



© 2011 RAS, MNRAS OOO.ITlfTOl 



Hadronic 7-ray images of SNRs. SN 1006. 



5 



3 GENERAL PROPERTIES OF HADRONIC 
IMAGES. NUMERICAL SIMULATIONS 

In order to narrow the parameter space, the images of SNRs 
in the present section are synthesized using the following 
set of parameters. Namely, we take 7 = 5/3, (5s = 0, s = 2, 
a = 1, K constant in time (i.e. b = 0), the energy of 7-rays 
e = ITeV. The actual representation for £'max(6o) used in 
our calculations is monoto nic function which corresponds to 
the time-limited model of Rey nolds! (|l998 l) with so called 
"gyrofactor" 77 = 1 (Fig. [1] [1 [5| or 77 = 5 (Fig.Elil); such 
form provides g « and Bp, max at the perpendicular shock 
to be 2 or 26 times the value at the parallel shock. The age 
of SNR and the number density are taken t — 1000 yrs and 
riHs = lcm~^ except of Fig. [5] where the densities are given 
in the captiorQ. Role of other choices of the above param- 
eters is visible from the approximate formula presented in 
Sect. H 

Radio and X-ray synchrotron images of SNRs de- 
pends directly on the distribution of magnetic field. Mag- 
netic field affects the TeV 7-ray images of SNR (due to 
inverse-Compton emission) indirectly, through modification 
of the downstream evolution of emitting relativistic elec- 
trons, namely, the larger the field the larger the radiative 
losses and, therefore, the thinner the region behind the 
shock occupied by the electrons able to emit 7-rays in the 
TeV band (Papers I and II). Magnetic field appears in the 
hadronic 7-rays only through dependence (if any) of the in- 
jection efficiency and the maximum energy on the shock 
obliquity. Therefore, the azimuthal variation of brightness 
in such images depends mainly on the two functions: /C(0o) 
and fmax(6o). 

3.1 Different models of injection 

Fig. [T] shows how Ok affects 7-ray image of SNR it has due 
to interactions of accelerated protons with thermal protons 
downstream of the shock. Different dependence of the in- 
jection efficiency on obliquity results in two different types 
of SNR morphology (Fig. [1]). Quasi-parallel injection cre- 
ates the 'polar-caps' structure (the lines of ambient mag- 
netic field cross the bright limbs) while the isotropic and 
the quasi-perpendicular injections are responsible for the 
'barrel-shaped' remnant (ambient magnetic field is paral- 
lel to limbs in this case). Increase of Ok results in decrease 
of the azimuthal width of the limbs due to the progressive 
luck of accelerated protons. 

The role of the aspect angle (po (between the ambient 
magnetic field and the line of sight) is shown on Fig. [2] 
Spherically-symmetric morphology transforms to bilateral 
one with increasing the angle if the proton injection depends 
on obliquity. 

3.2 Dependence on the maximum energy 

The ratio of -Emax at the perpendicular and parallel shock is 
2 on previous figures. Let the maximum energy of protons 
increase more rapidly from parallel to perpendicular shock, 

^ parameter really influencing a 7-ray image of SNR is a prod- 
uct of t and rtHsi Eq. 1 II6II . 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 



r/R r/R 

Figure 3. Role of -Emax in hadronic 7-ray images of Sedov SNR. 
£^max|| = 10 TeV (left, = 0.82), EmaxH = 1000 TeV (right, 
Em = 0.0082). Isotropic injection, £max(7i"/2)/£max(0) = 26, o = 
1, aspect angle <j>o = 90°. 



in particular in 26 times as on Fig. [S] In this case, the az- 
imuthal variation of i?max is in general more prominent. 

The thickness of limbs depends on two factors. First, 
on the ratio em = -Epm/-Ep,max|| where -Epm is the energy 
of proton which effectively gives the most contribution to 
7-rays with energy e; it is given by Eq. (|D5P : once we are 
interested in e > 1 GeV, it may be simplified to 

£pm(£) =eC/K = 8.2e. (18) 

where we used f = 1.4 (Appendix |D} . Namely, if em ^ 1 
then the obliquity variation of the maximum energy does not 
affect the image of SNR (Fig. [Sjl . Following this conclusion, 
we would like to note that the role of -Emax may be promi- 
nent in the analysis of the Fermi observations (e ~ 3 GeV) 
only if -Ep,max is smaller than ~ 25 GeV that is unreason- 
able. Therefore, the role of -Ep,max may be neglected during 
analysis of observations in GeV 7-rays. 

The second factor which affect the azimuthal thickness 
of limbs is how quickly the proton energy spectrum N{Ep) 
ends (given by the parameter a in our model). Fig. Udemon- 
strates this effect: the larger a the more rapid decrease of 
the high-energy end of the proton spectrum and the smaller 
the azimuthal extension of the limbs. 



3.3 Role of density of the target protons 

Energy losses of relativistic protons due to collisions with 
thermal protons are proportional to the density of the tar- 
get protons, Eq. (jlOp . They were negligible in the calcu- 
lations presented above because the number density was 
whs = lcm~^. The coUisional losses of protons are impor- 
tant for much higher densities as it is evident from Eq. H16p . 
The role of the pre-shock density in hadronic 7-ray images 
of SNRs is similar to the role of the magnetic field strength 
in the hard X-ray images. Namely, the larger the density 
the higher the coUisional losses of accelerated protons and, 
therefore, the thinner the radial profiles of 7-ray brightness 
(Fig.O. This effect may be used for estimation of the target- 
proton density from thickness of 7-ray rims, once the 7-ray 
observations reach necessary resolution. 
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Figure 4. The same as on Fig.|3]for a = 0.5 (left) and 2 (right); 



4 APPROXIMATE ANALYTICAL FORMULA 
FOR HADRONIC 7-RAY IMAGES 

The approximate formulae (valid close to the shock front) 
for the radial and azimuthal profiles of the surface bright- 
ness of Sedov SNR due to leptonic emission (synchrotron 
radio, X-rays and inverse-Compton 7-rays) are presented in 
Papers I and II. The same approach is used here to derive an 
approximation for the surface brightness profiles of hadronic 
7-rays arising from internal structures inside adiabatic SNR. 
Namely (Appendix ID]). 



Sp{ip, q) oc ^(9o,efi) exp 



/pp(eo,ofi,e)(i9) 



where 



(20) 



1 - 



pjCT(Kpp+2) 



1 - ^<^(«PP + 1) Kpp +2 



Co-p 



: 3fe/2 + (5 + s)Kad + 1/0- - f + Ccrpp(s - 1), V = Kad + 

- 3g/2, Kad ~ 1, (T the shock compression ratio, p = 



p/R- Note, that the azimuthal variation of Ipp arises only 
from obhquity dependence fmax(6o,cff )• 

The effective obliquity angle So, off, the azimuth tp (mea- 
sured from the direction of ISMF in the plane of the sky) 
and the aspect angle (j>o are related as 



cos 0o,cff (<y5i 't'o) ~ COS if sin <j!>o 



(21) 



Approximation p9p is compared with numerical calcu- 
lations in Appendix [D] The formula are rather accurate in 
description of the brightness distribution close to the shock. 
It restores all the properties of the surface brightness pro- 
files of Sedov SNR, revealed in the numerical simulations, 
including dependence on the aspect angle. They are quite 
useful for qualitative analysis of how different factors af- 
fect the hadronic image. They may also be used as a simple 
quantitative diagnostic tools for hadronic maps of SNRs. 

For example, the case when the number density is not 
high, C, <^ 1. Then, for known a, the thickness of the 7- 
ray rim depends mostly on the compression ratio a and 6; 
also on Q, Eni and q; the later could be taken zero for an 
assumption that i5max is limited by the time of acceleration. 
An opposite case, the product nH(t/1000yrs) > 10^. Then, 
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Figure 5. Radial profiles of 7-ray surface brightness due to 
hadronic emission for different densities njjs of the target protons: 
lcm-3 (line 1), lO'^cm-^ (line 2), 10^ cm'^ (line 3), lO'^ cm'^ 



(line 4). Isotropic injection, £max(8o) = const, 



1000 



TeV, 0o = 90°, azimuth angle (measured from direction of the 
ambient magnetic field) 1/3 = 0°. 5max is the peak value of Sp{p). 



of the 7-ray rim may be used to estimate the post-shock 
number density. 

It is merit to note that our approximate formula re- 
flects also some effects of the non-linear acceleration theory. 
Namely, in case of efficient proton acceleration, a prominent 
fraction of the kinetic energy of the shock goes to the rela- 
tivistic particles. If so, plasma may be described by the adia- 
batic index 7 smaller than 5/3. In our formula, the adiabatic 
index appears through the compression ratio a (Eq. IBT9)| 
and the parameter Kad (Eq. IC3I though Kad is close to unity 
in the range 7 = 1.1 -;- 5/3). Features on images are radially 
thinner for smaller 7. Eventual amplification of magnetic 
field would be accounted through the larger magnetic field 
compression ratio; however, this ratio does not appear in 
our formulae because the hadronic images of SNRs do not 
depend on the strength of magnetic field. 



5 DISCUSSION AND CONCLUSIONS 

The spectrum of the TeV 7-rays from SN 1006 may be inter - 
preted as leptonic as hadronic in origin ([Acero et al.ll2010l ). 
The data in the GeV 7-ray range is expected to constrain 
further this ambiguity. The pattern of surface brightness of 
SNR also contains important information. In the present pa- 
per, we consider possibility that some 7-ray images of SNRs 
may be hadronic in origin. Namely, we are interested in the 
images the Sedov SNRs would have if accelerated protons 
interact with the thermal protons downstream of the shock. 
The model to synthesize maps of the 7-ray surface bright- 
ness of adiabatic SNRs in uniform ISM and uniform ISMF is 
developed. It includes parameterized description of surface 
variations of parameters characterizing the injection and ac- 
celeration of protons as well as evolution of relativistic pro- 
tons downstream. The later considers both adiabatic losses 
of energy and losses due to inelastic collisions. Collisional 
losses are non-negligible only in case of the large number 
density of target protons. For example, energy losses of pro- 
ton with energy 30 TeV due to pp-interactions are compara- 
ble to the radiative losses of electrons with the same energy 
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in magnetic field 30 fiG, is the number density of the target 
protons is ~ 10^ cm"^. In case the shock moves in the ambient 
medium with smaller number density, one can consider only 
the adiabatic energy losses of accelerated protons to model 
an image and volume-integrated spectrum of SNR. 

The radial thickness of the hadronic 7-ray rim in case 
of SNR interaction with the large-density ambient material 
may be used to estimate the density. 

The azimuthal variations of the surface brightness of the 
shell-like adiabatic SNR in hadronic 7-rays is a consequence 
of the obliquity dependence of the injection efficiency and 
the maximum energy of accelerated protons. The orientation 
of the pre-shock magnetic field is the only way the magnetic 
field infiuences hadronic images of SNR. If magnetic field is 
highly turbulent everywhere before the shock of SNR evolv- 
ing in the uniform ISM then SNR should look as a ring. 
Nonuniform ISMF (different strength over the SNR surface) 
does not able to provide any deviation from the ring pattern 
if information about orientation of ISMF is lost due to turbu- 
lence. In contrast, ISM with nonuniform density distribution 
might provide patterns other than ring, since K (x n. How- 
ever, it is unlikely that quite symmetrical bilateral structure 
(cis in SN 1006) may appear due to nonuniformity of ISM 
density because the structure of ISM should be quite special 
(like a tube, which is not observed around SN 1006). 

Bilateral pattern in hadronic image of the shell- like SNR 
may therefore be considered as a sign of an ordered ambient 
magnetic field. If so, the limbs are due to azimuthal variation 
of either /C or £'max. The most contribution to 7-rays with 
e ~ 1 TeV gives protons with energy _Ep ~ 8 TeV. There- 
fore, if -Bp, max ^ 8 TeVQ then the bilateral pattern of the 
hadronic TeV image of SNR reveals the variation of the in- 
jection efficiency of protons; in fact, Eq. (|19|l simplifies to 
Sp{^p) oc ?(0o,cff (Vi "^o)). In this case, the hadronic GeV im- 
age has to be the same as TeV 7-ray image (differences in 
GeV and TeV 7-ray images of SNR are signs of the smaller 
Sp,max or the contribution from the leptonic emission) and 
the azimuthal variation of the hadronic 7-ray brightness may 
be used to derive the obliquity dependence of the proton in- 
jection efficiency. Unfortunately, large errors in the present 
7-ray data on SN 1006 prevent us from possibility of such 
an estimate. 
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APPENDIX A: ENERGY LOSSES OF 
PROTONS DUE TO PION PRODUCTION 

In order to obtain the energy losses of a "single" proton due 
to pion production, one needs to integrate the pion source 
sp ectra over energie s of pions Et^. The pion source spectra 
is (|Schhckeiserll2002l ) 

p[E^,E^) = cnnE^^^^^^%^H{E^-E,^), (Al) 

where c, nn are speed of light and the number density of tar- 
get protons, dapp{ETy, Ep)/dETr the differential cross-section 
for the interaction of two protons, Ep the energy of incident 
proton, H denotes the Heaviside step function and Eth the 
threshold energy of interaction. The energy losses due to 
inelastic proton collisions is 



( dEp 
\ dt 



dE^p{Ep,E^) 



(A2) 



where the factor 3 accounts for the production of and 
n~ mesons respectively. In the (5-func tion approxiraation, 
the differential cross section is given by (|Dermeilll986al : lMoril 
1 19971 : lAharonian fc Atovan| [2000") 

dcjpp{E-^, Ep) 



dE„ 



= app{Ep)S{ETr - K-Bp,kin), 



The proton collision losses are 
_ (dEp 



V dt 



= 3KCnHCrpp(i5p)i?p,kin 



(A3) 



(A4) 



where k = 0.17 (|Aharonian fc Atovanllioo d) and therefore 
3k = 0.51 tha t restores the coe fficient of inelasticity / « 0.5 
(Eq. (3.13) in lAharonia"nll2004l '). 



APPENDIX B: EVOLUTION OF THE PROTON 
ENERGY SPECTRUM DOWNSTREAM OF 
THE SHOCK IN SEDOV SNR 

We assume that relativistic protons are confined in the fiuid 
elem ents which advects them from the region of acceleration 
iRevnolda .1998. 1^. An individual proton loses energy due 

^ This means that we do not account for diffusion in the present 
consideration. This is the reasonable assumption for protons with 
Ep ~ 1 TeV (and smaller energies) which give most contribution 
to emission at TeV (GeV) photons respectively, and for the case 
when the lengthscale of diffusion is proportional to gyroradius. 
In fact, the gyroradius of the relativistic proton is the same as 
gyroradius for electron with the same energy; the lengthscale of 
diffusion for electrons with energy Ee ~ 1 TeV (and even for 
higher ene rgies) are sim ilar or smaller than the lengthscale for 
advection jBallet ll2o6^) . 



to inelastic collisions and adiabatic expansion. Downstream 
distribution of the target protons is pro portional to the den- 
sity distribution given bv ISedovl (| 19591 ) solution, namely it 
is p(r) — p{r)/ps where index "s" marks the immediately 
post-shock value. The pion production losses are given by 



Ep 



3Kcnji^sp{r)app{Ep)Ep 



(Bl) 



(here, we neglect difference between total and kinetic ener- 
gies of proton because we are mostly int erested in Ep ^ 
IGeV). The adiabatic losses are given by (iRevnoldjliggj ) 



_• _lEpdV__Epdl 

^P.ad - 3 ^ - 3^ rf^ • 

The equation for losses is 



(B2) 



dEp 



+ 



Ep dp 



3K,cnHBp{r)app{Ep)Ep. 



dt 3p dt 

In therms of w = Ep/p^^^ (| Reynolds! 1 1991 ). Eq. ((BS]) is 



(B3) 



— = -cxpwapp{wp ), 



(B4) 



where Ci = 3 x 10 ^^Kcnns cgs. 

Let us take the cross section Opp in the simple form 
jAharonian fc Atovanll2000l ) 



'^Pp(-Ep) = 28.5 -h 1.8 In (ca-Ep) mb 

where ca = (lGeV)~\ 

The solution of this equation is 



In 



c-iEp\ _ ln(c2£pi) + In/ 



51/3 



(B5) 



(B6) 



where Ep-^ is the initial proton energy produced on the shock 
at time ti. 



lit) — exp 



p{t) — exp 



ci / p{t")dt" 



-p (28.5 -h 1.8 In p 



l/3\ 



Ci / q{t')^l{t')dt' , (B7) 

(B8) 

(B9) 
(BIO) 

It follows from (|B6P that Epi is related to proton energy Ep 
at the time of interest t as Ep — Epi£!^^£pp where 

Cad = P 

represents the adiabatic losses and 
£pp = (ca^p)'""/ 

represents the energy losses due to inelastic collisions. 

Self-similarity of the task allows us to write p and / 
in terms of the normalized Lagrangian coordinate a = a/R 
rather than in terms of time t. Namely, I{d) and p ( a) are 
dimensionless functions and, with the use of lSedovl (|1959D 
solutions for uniform ambient medium, dt' /t = {5/2)x^^^dx 
where x{t') = R{t')/ R{t), and therefore 



(Bll) 



(B12) 



p(a) 



exp 



3/2 
X ' p 



dx 



(B13) 
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I{a) = exp 



dx 



(B14) 



where C = 5tei/2 = 1.21 x 10"''t3nHs, U = t/lOOOyrs. It is 
clear from here that fpp is effective only where the number 
density of target protons is large, at least tihs ~ 10® cm~'^. 
Thus, relations between energies and energy intervals 

are 



F 



dEpi 



fi{a) 



dEp £^i{dY^-^)£pp{E,dy 



(B15) 
(B16) 



Let us assume that, at time U, a proton distribution has 
been produced at the shock 



Ep 



where a is constant. The conservation equation 
a^dadEr^i 



N{Ep,a,t) = N{Epi,a,ti) 
where 

^_ ps 7 + 1 



ar^drdEp ' 



7 - 



(B17) 



(B18) 



(B19) 



is the shock compression ratio (index "o" marks the pre- 
shock value), and the continuity equation po{a)a?da — 
p{a,t)r^dr yield that the energy spectrum evolves down- 
stream as 



Np{Ep,a,t) ^Kia^tJEp" fi{a)£pp{Ep,dy 



X exp 



iJpa3'/2 



Ep 



>=£ad(a)'^(")£pp(Sp,a) 



(B20) 



with K{a,t) = A'a(fi)pfad(o)'''"'*''"^^ the evolution of K is 



self-similar downstream for Kb (x V 



K{d) = K{a,t)/KB{t) = d^''^^p{d) 



0+M{a){3-l)/3 



(B21) 



We note that if the cross section app ~ const then jJ. — 1 
and q = o-ppp. 



APPENDIX C: APPROXIMATIONS FOR 
EVOLUTION OF SOME PARAMETERS 
BEHIND THE SHOCK 

We interested in approximations of some function x down- 
stream and close to the shock, in the form 



X{a) ^ a". 
By definition 

a dx{a) 



X{a) da 





91nx(a) 




9 In a 



(CI) 



(C2) 



This approach yields for adiabatic losses (jPetruk et al.l 
l2011bl ) 



57+13 
3(7 + 1)^ 



(C3) 



it is valid for r > 0.8 with error less than few per cent. 
The value of Kad = 1 for 7 = 5/3 and is close to unity for 

7 = 1.1-5/3. 

Applying ([C2| to Eqs. (|BT3|) and (|BT4)) one obtains: 

^i(a) a-'-*«, I{d)^a^''-^'^. (C4) 



With these approximations, Eq. (|B12p are approximately 

£:pp(£p,a)^a'^"pp(-^p' (C5) 

where the cross section is in units of mb and we used the 
property s^"" ~ a^'"'^^ which is valid for d ~ 1. Approx- 
imation (|C5|) is independent of the adiabatic index 7. 



APPENDIX D: APPROXIMATE FORMULA 
FOR PROFILES OF THE HADRONIC SURFACE 
BRIGHTNESS OF SEDOV SNR 

Dl Derivation of the formula 

1. The 7-ray emissivity due to meson decays is 



1. ine 7-ray emissivity 
(|Aharonian fc Atovanll200oi ') 

f qTr(ETr)dE-^ photons 
'^'^'>-' J ^Ei - mic^ cm^seV 

the pion emissivity 

g^iJ,) = £!^app {E'p)N{E'p). 



(Dl) 



(D2) 



On average, energy transferred from protons to pions is 
Ejy = fciJp^kin, where k = 0.17 accounts also for contribu- 
tion from ry-mesons; therefore an average energy of protons 
mostly responsible for creation of pions with energy i?^ is 

Ep = rripc^ + E^/k. (D3) 

The minimum energy of pion to create photon is 

£^,min = e + mlc^/Ae. (D4) 

Substitution (|D4|) into (|D3P results in expression for E'p -^^-^, 
an average minimum energy of proton to create photon with 
energy e. 

Let us introduce an effective energy of protons 
£'pin(e) = C^'p.min(£) which gives most contribution to 7- 
rays with energy e: 



f / rn^c 2 

Epraie) = - e H h unipC 

^ ^ ' k\ 4e 

and re-write Eq. (|D1|I in the form 



q^{e) = '^^app{Ep^)N{Ep^)Fp{e,i) 
where 



F{e,i) 



E' . (e) 



(D5) 



(D6) 



(D7) 



F(e,5) ^ 1 is provided by the values ^ = 1.42, 1.40, 1.44 for 
s = 2, 1.8, 2.2 respectively, in a wide range of Ep (Fig. IDlj ). 
To the end, Eq. (|D6|) becomes 



(D8) 



g^(e) ~ crpp(_£;pni)A''(-Epm) 
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s=2.0;a=1 .0:4=1.42 
s=1.8;a=1. 0:^=1.40 
s=2.2;a=1 .0:4=1.44 
s=2.0;a=0.5:4=1.42 
s=2.0;a=2.0:4=1.42 



-5 -4 



-2 -1 

Ig E.TeV 



Figure Dl. Function F(e), Eq. | |D7|1 . for few sets of parameters. 
Ep,max = 10^^ eV. Rapid deviation the function from unity at 
energies above e ~ 0.1-Ep,max corresponds to the rapid decrease 
of 7-ray emissivity. 



We assume that, close to the shock front, the cross section 
iTpp ~ const. The distribution of target thermal protons is 
proportional to the plasma density, nn oc n(a). Therefore, 
Eq. (|D8|) may approximately be represented with 

q-,{e)<xn{a)N{Ep^,a). (D9) 

2. The energy of protons in a given fluid element at 
present time was _Epi at the time this element was shocked; 
they are related with HB15|I . The behavior of the factors fad 



and £pp close to the shock is given by Eqs. HC3|I and (|C5 
The distribution N(Ep) is therefore approximately 



N{Ep,eo) oi]C{Qo)K{a)E-' 



Epa-1' 



Ep^raax, ||^max 



(eo) 



(DIO) 



1), 



(Dll) 



X exp 

where V'2(-Bp) = (app{s - 
= "^ad + C'^pp - 3g/2 

and we use fJ, ~ 1 for app ~ const. 

3. An 'effective' obliquity angle 0o,c£F is related to the 
azimuthal angle ip and the aspect angle (/)o as (Papers I and 
11) 



cos Oo,cff (</5, 0o) = cos ip sin ( 



(D12) 



where the azimuth is measured from the direction of ISMF 
in the plane of the sky. 

The surface brightness of SNR is 

S{g,p) = 2[ q,{a) "J^^ , (D13) 



where g is distance from the centre of the SNR projection, 
fa is the derivative of f(a) in respect to a. With the use of 
Eqs. (jD9|) and (|D10|l . the variation of the meson-decay 7-ray 
brightness is approximately 



Sp oc ?(9o,cfi) exp 



.(e) 



^-max ( o , cff ) 



/pp(eo,eff,^) (D14) 



where em = -Bpm/-Ep^max,| 



I a(e) 



■ exp 



-(a 



(D15) 



4. Let us approximate /pp. Close to the shock front, the 
density distribution is 



-3k 

n ~ a 



and Eq. HB21|) yields therefore 

In addition l|Petruk et aLlbOllbT l: fa ^ (l/cr)a 
exp {—qt,{a~°"^ — 1)) 1 — qtaip{l — a) and 

f 1 



1/(T-1 



The integral of interest is therefore 



(D16) 
(D17) 

(D18) 



T{ftpp + 1) 



e° qV- / 1 - e'^(''PP + 2) Kp 



(Dig) 



1 — gO-(Kpp+l) 



where 



36 



+ (5 -I- s)Kad + 1 + C'^pp(s - 1). 



(D20) 



— PP o • - ^ 

The formula Eq. (fPTI)) gives us the possibility to ap- 
proximate both the azimuthal and the radial brightness pro- 
file for Q close to unity. 

D2 Accuracy of the approximation 



Fig. ID2I demonstrates accuracy of the approximation HD14|) 
(upper and lower panels for azimuthal, radial proflles). Our 
calculations shows that this approximation may be used, 
with errors less than ~ 20 -f- 30%, in the range of q from 
0.9 to 1, for azimuth Lp where e/-Ep,max < 0.5, 50, 0.2 for 
Q = 1, 0.5, 2 respectively (Fig. IDl]) . Accuracy of approxi- 
mation is higher for smaller 0k (because the injection term 
dominates in the azimuthal variation of the brightness) and 
for smaller variations of £max (the role of .Smax becomes how- 
ever less prominent for smaller ratio e/iSmax)- 
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0.9 0.92 0.94 0.96 0.98 1 0.9 0.92 0.94 0.96 0.98 1 0.9 0.92 0.94 0.96 0.98 1 

p/R p/R p/R 

Figure D2. Azimuthal (upper panels) and radial (lower panels) profiles of the hadronic 7-ray surface brightness Sp (solid lines) and 
their approximations IID14t (dashed lines). Calculations are done for (fiQ = 90°, quasi-parallel injection with 0k = 7r/2, other parameters 
are the same as in Sect. |3] Models of i?niax: blue line - fmax = const; black line - the smooth mo notonic variatio n which provide the 
variation of ^max in 3.3 times with azimuth (actual representation used is the time-limited model of iRevnoldsl lll998l ) with rj = 1.5). The 
ratio of the photon energy to the maximum proton energy is £/^/p,max — 0-1 (l^ft); ^/^max — 0.01 (middle), s/ii/max — 0.001 (right). 
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